USTC-ICTS-10-02 



I. INTRODUCTION 



Cosmic age test in inhomogeneous cosmological models mimicking ACDM on the light 

cone 

Mi-Xiang LaV' 2 'Q Mia o Li, 3 - 2 'Q Xiao-Dong Li, 1 ' 4 ' 2 '© and Shuang Wang 1 - 2 '! 
1 Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China 
2 Key Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, China 
s Kavli Institute for Theoretical Physics China, Chinese Academy of Sciences, Beijing 100190, China 
^Interdisciplinary Center for Theoretical Study, University of Science and Technology of China, Hefei 230026, China 

The possibility of reconstructing a spherically symmetric inhomogeneous Lemaitre-Tolman-Bondi 
(LTB) model with ACDM observations has drawn much attention. Recently, an inhomogeneous LTB 
model having the same luminosity-distance and light-cone mass density of the homogeneous ACDM 
model was reconstructed. From the Wilkinson microwave anisotropy probe 7-year measurements 
together with other cosmological observations, we calculate the cosmic age at our position in this 
LTB model, and obtain a constraint ti,TB < 11.7Gyr at la confidence level. We find that this 
result is, although 2Gyr younger compared with the age of the homogeneous ACDM model, still 
(N ■ within la agreement with the constraint of cosmic age given by current astronomical measurements. 

We expect that in the future with the help of more advanced observations we can distinguish the 
reconstructed inhomogeneous LTB model from the homogeneous ACDM model. 

u 

Or 

6 

<*— » , The problem of dark energy has become one of the most important issues of the modern cosmology since the 
observations of type la supernovae (SNe la) [l[ first indicated that the universe is undergoing an accelerated expansion 

i— —i \ at the present stage (if assuming that the universe is described by the Friedmann-Lemaitre- Robertson- Walker (FLRW) 
model). Many cosmologists believe that the identity of dark energy is the cosmological constant which fits the 
observational data very well. However, one still has reasons to dislike the cosmological constant since it suffers 
from the theoretical problems such as the "fine-tuning" and the "cosmic coincidence" puzzles [2j]. Thus, a variety of 
proposals for dynamic dark energy have emerged. For example, the "scalar field" model Q has been explored for a 
, long time. Besides, the "holographic dark energy" models [J, which arise from the holographic principle of quantum 

C^) ' gravity theory, has also attracted much attention. 

, There are also some other theoretical approaches to explain the current cosmic acceleration. For example, it is 
argued that the acceleration of the universe may signify the breakdown of Einstein's theory of general relativity § ■ 
Another interesting idea [f| is based on the assumption that the universe is described by the spherically symmetric, 

t— I . inhomogeneous Lemaitre-Tolman-Bondi (LTB) metric [1] @- Recently, some authors further proposed a possibil- 
ity of mimicking the cosmological constant in an inhomogeneous universe (Toj [ill ] (l2j [L3| . The idea is that, since 

• i-H , cosmological observations are limited on the light cone, it is possible to reconstruct an inhomogeneous cosmological 
model (indistinguishable from the homogeneous ACDM model) to explain the cosmic acceleration without a cosmo- 
logical constant. The formalism of reconstruction was developed in [1(| and was applied to the ACDM model in 
[ljj . In [l2| . the authors constructed a spherically symmetric, inhomogeneous cosmological model reproducing the 
luminosity-distance and the light-cone mass density of ACDM model up to z = 2. 

In this work, we focus on the reconstructed inhomogeneous cosmological model and investigate whether it is con- 
sistent with cosmological observations. Since in this model ACDM observations are exactly reconstructed on the 
light-cone, we should seek for some observational test not limited on the light-cone to distinguish it from the standard 
ACDM model. 

Fortunately, we find that the age of the universe is an appropriate touchstone. The cosmic age, which depends on 
the evolution of the universe at a comoving position, contains information not limited on the light-cone (the age is 
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uniform in the ACDM model but may be dependent on the position in an inhomogeneous cosmological model) . Thus 
it may reveal the discrepancies between the ACDM model and the reconstructed inhomogeneous model. On the other 
hand, it is rather convenient to use the cosmic age to test the validity of a specific cosmological model. To do this 
one may just compare the result with the age of some old objects in our universe For example, to be consistent 
the age of our universe in our position must not be younger than the age of the oldest stellar in the Milky Way. 

This paper is organized as follows. In Sec. HI following the procedure of [13], we introduce the inhomogeneous LTB 
model and explain how to reconstruct ACDM observations in this model. In Sec. IIII1 we calculate the cosmic age at 
our position and compare the result with the age of some old objects. We summarize in Sec. IIVI 



II. LEMAITRE-TOLMAN-BONDI MODELS 

In this section, following the procedure of [l2j , we explain how to reconstruct an inhomogeneous LTB model having 
the same luminosity distance and light-cone mass density of the homogeneous ACDM model. 

The LTB models are spherically symmetric cosmological solutions to the Einstein equations with a dust stress-energy 
tensor. The general metric for the LTB models in a synchronous comoving coordinate takes the form, 

ds 2 = -dt 2 + R ' 2{ ^\ dr 2 + R 2 (r, t)dtf. (1) 
1 + j3{r) v ' w 

Following [l(| H3], we use the prime superscript to denote d/dr, and the overdot to denote d/dt. Noticing that 
here r is a dimensionless coordinate, while R(r, t) has the dimension of length. The Robertson- Walker metric can 
be recovered after performing R(r,t) — > a(t)r and f3(r) — > — kr 2 . Solving the Einstein Equations one obtains the 
generalized "Friedmann Equations" for R(r,t) and p(r,t), 



np(r, t) 



a'{r) 
R 2 (r,t)R'{r,t)' 



(3) 



And the photon radial null geodesic equation for t(r) is found directly from the LTB metric (we are only interested 
in the past light cone), 



dt(r) 



R'(r,t(r)) 



For convenence we denote quantities on the light cone by a hat. Thus we have, 

R(r, i») = R; R'(r, t») = R'; p{r, t(r)) = p. 



(4) 



(5) 



Then let us focus on the reconstruction procedures. The method was discussed by Mustapha, Hellaby, and Ellis in 
1997 [l(J, and was recently applied to the ACDM model in 11] pjj (A related formalism was developed by 15] [l6j). 
Following their procedure, we take advantage of a coordinate freedom and rescale r so that on the light cone, 



# = fT o -yi + 0(r). 



(6) 



The corresponding coordinate transformation is dr\ = H( ' dR ^ dr dr. In the following, for simplicity we will use r to 

\/l+/3(i-) 

denote r\. The redshift takes the form (see Sec. 2.3 in pjj]) 



dz(r) 
dr 



(l + z) 



R' 



Vi + W) 



(J) 



Following the procedure of 12|, we reconstruct dh{z) and p(z) of the ACDM model. To do this we require that on 
the light cone, 



(l + z) 2 R(z) = (1 + *) 



dz\ 



o -Hacdm(zi 



-, p{z)dVhTB = PM.ACDMrfVACDM, 



(8) 
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where Pm,ACT>m(z) stands for mass density in the ACDM model, and -Hacdm {z) stands for the Hubble constant in 
the ACDM model. These quantities take the forms, 



Pacdm(z) = 3r2 TO Mpi7 ACDM (z) 2 , H ACBM = H oy /n m (l + z) 3 + Q. A . (9) 
For simplicity we only consider a flat ACDM model. We use the notations, 

On = ^4nT, ^A = ^7nT, Pc = 3M p 2 ff 2 , (l m + fi A = 1. (10) 



Then from Eq. ([5]) it is straightforward to derive the following expressions 

1 f z dz\ 

l + zj i?ACDM(2l) 

H- 1 R 2 {z)np(z) dr = 3flmH v [ r rfzi 1 2 (12) 
dz Hacbm(z) Jo Hacbu(zi) 

Furthermore, the following three equations can be obtained by solving the corresponding z, a(r) and /3(r) (Eqs. 
(19-21) in dl) 

? = (l+*)^f^, (13) 
da 1 , , , . r 1 / a x dR 1 

It should be emphasized that for these values of u(z) and (3(z), the corresponding LTB model exactly reproduces 
the luminosity-distance relation and light-cone mass density of the ACDM model. For convenience, in the following 
context we will use "LTB-ACDM model" to denote this reconstructed LTB model. 

It should be stressed that since the LTB-ACDM model is reconstructed by mimicking ACDM model, it has the 
same luminosity-distance-redshift relation and light-cone mass-density-redshift relation as the homogeneous ACDM 
model. Thus when estimating some physical quantities of the LTB-ACDM model from cosmological measurements 
one can just use the obtained values of f2 m and Hq from the fit of the ACDM model (i.e. it is not necessary to impose 
a special constraint on the LTB-ACDM model). For example, in [l2| the authors just take fi m =0.3 to mimicking a 
flat ACDM with mass ratio p m o/PcO — 0.3. 



III. COSMIC AGE AT OUR POSITION 



In the first subsection, we derive the expression of the cosmic age in the LTB-ACDM model. Then in the sec- 
ond subsection, we estimate the age from cosmological observations, and compare it with that in the homogeneous 
ACDM model. Finally, in the third subsection, we discuss the validity of the LTB-ACDM model by considering the 
astronomical measurements of the age of some old objects in our universe. 



A. Cosmic Age in the LTB-ACDM 



First of all, we investigate the properties of Eqs.([2^).(l3^). (|lHI'T5l) at r = (corresponds to our position in the universe). 
Since the function R(r, t) represents the diameter distance, one expects R{r, t) — > when r approaches zero. So R{r, t) 
can be expanded near r — as, 



R(r,t) = R 1 {t)r + 



(16) 
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Then from Eq. (| 14[) we obtain 



a(r) = a 3 r d + 0{r) = for 4 + Rx(t) = J fa + (17) 

From Eqs. © (HH) (HU) , we expand i?ACDM(z), -R(z) to the leading order and substitute them into Eq. ([H]) to 

obtain 03 

Hacdm(z) = -ffo + -, r = z + R(z) = -^- + a 3 = ^ m H^ . (18) 
The calculation of /?2 is also straightforward. From Eqs. (fT4"))(|15p it follows that 

. . ,1 r 1 / Ol \ (1R 1 X 2 , 

Substituting fl'(z) = Hq 1 + R[z + R' 2 z 2 + ... into Eq. we obtain 

P(z) = {HqR'i - n m )z 2 + ... = (HqR'i - a m y + ... (20) 

Notice that R' 2 does not appear in the expression, so we just have to expand R! to the first order. Using Eqs. (fTT|) (|T3 )) . 
it follows that, 

dR = dRdz = _ 1 H ACnM (z) r dzx J_ =H i_ H i z+ ^ R [= H^. (21) 

dr dzdr 1 + 2 H J H A cnu(zi) H ' 10 V ; 

Combining with Eq. (|20l) we obtain 

P2 = l- On- (22) 

Next we calculate the age of the universe (to — tss) at r = 0. From Eq. (IT71) it follows that, 

/•to /-Hi (to) / t \ 

to-t BB = dt = -^Aii, (23) 

JtBB JRl(tBB) R\{t) 

where the upper and lower bounds of the integral can be determined by Eq. ([3]), 

Finally, combining the above with Eqs. ([5)1 (ITS)) , we obtain the age of the universe at r — in the LTB-ACDM model, 

f"" 1 dRi , , 

to-<ss= / , ==■ (25) 

SI a + il m H I R\ 



n 



B. Estimate tLTB from Cosmological Observations 

Next we estimate the cosmic age in the LTB-ACDM model and compare the result with that of the ACDM model. 
These two models are both determined by two parameters, the present matter ratio Q m , and the Hubble constant Ho- 
For convenience, let us denote the cosmic age in the LTB-ACDM model and ACDM model at r — as iLTs(^m, Ho) 
and tACDM^m, Ho) ■ We have [we perform a parameter transformation r — H^ 1 /(I + z) in £acdm] 

f H ° X dRi 

t LTB (n m ,H )= . = p==== , (26) 



^^/(HoR^ + l-nJ 

f +ca dz _ f 11 " 1 dr 

tACDM(ll m ,Ho) — / 5~TT7i — V~ / /„ 77= , 7= = ?t== = ■ I 2 ') 

Jo H A (z)(l + z) Jo y/Sl m /(H r) + (1 - n m )(H r) 2 
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The only difference between £acdm and tLTB lies in the second term in the square root. Since H Q r < 1 (the upper 
bound of the integral is H^ 1 ), it is clearly that tLTB is smaller than ^acdm given the same set of parameters. 

Notice that tLTB^m, Ho) is proportional to 1/Hq. To prove this one can perform a coordinate transformation 
i?i = HqRi, which yields that, 

t LTB (n m , Ho) = i- I = =. (28) 

So we have, 

9 -t LTB {n m ,H ) = -^-t LTB (n m ,H a ) < 0, (29) 



H ° Jo Jn m /Rl + 1 - fir, 



v ' ' Ho 

9 -t LTB (n m ,H ) = -\ f H ° ( ^ + 1 - fi ro )- 3/2 (-L - 1) dr < 0, (30) 



dfl m 2 Jo 

Clearly, tLTBi^m, Ho) has larger values with smaller values of fl m and Hq parameters. One can easily verify that the 
conclusion is the same for t A CDM(fim,-Ho)- 

Now we are ready to calculate the specific values of tLTB with given values of Ho and VL m . Here we refer to the 
result of the seven-year Wilkinson microwave anisotropy probe (WMAP) observations [l7[. From the constraint from 
"WMAP7+BAO+# " [13 [III El the WMAP collaboration provides the best-fit values of fi A and H together with 
their la uncertainties, 

fl A = 0.728±°;°^, H = 70.4±ij km/s/Mpc. (31) 

From their result we can put a constraint on tLTB and £acdm at la confidence level (CL). The result is, 

t LTB = 11.4 ± 0.3 Gyr, t ACD M = 13.8 ± 0.5 Gyr. (32) 

It is found that tLTB is about 2Gyr younger than £acdm- At the la CL we obtain the upper limit of the age of the LTB- 
ACDM model thTB < 11.7Gyr, with the set of the smallest values of parameters fi m = 0.257 and Ho = 69.0km/s/Mpc 
(In fact this result is overestimated since we ignore the degeneracy between VL m and Ho- The upper limit value of 
^ltb should be smaller, or at least as small as 11.7Gyr). 

It should be stressed that to be strict the previous estimation oHltb is not appropriate, since we have assumed that 
the WMAP and baryon acoustic oscillations data could be used to constrain the LTB models. In fact, the issues of 
CMB and structure formation in the LTB scenario are rather complicated and have not been clearly investigated. To 
avoid this problem one can put constraints to f2 m and Hq from SNIa observations. The recent observations of SDSS- 
II (Sloan Digital Sky Survey II) [20j] and Hubble Space Telescoped!] show that fl m > 0.224, H > 70.6km/s/Mpc 
in la CL. From their result we find a upper limit tLTB < 11.6Gyr, which is the similar with our previous result 
t LTB < llJGyr. 



C. Discussions of the Validity of the LTB-ACDM Model 



The low limit to the cosmic age can be directly obtained from estimating the age of some old objects in our universe 
21] [22] 23}. As an example, based on white dwarf cooling the authors of [2^] get a result of 12.7±0.7Gyr. Compared 
with the result of the previous subsection t^TB < 11.7Gyr it seems that the LTB-ACDM model is inconsistent with 
their measurements. However, one should not conclude hastily. The reason is that the result of [22| is subject to 
larger uncertainty. The uncertainty due to calculations of the white dwarf cooling is difficult to estimate, and in their 
result corresponding errors are not included. In fact, the authors of [22[ argued that systematic uncertainties are 
likely to be at least as large as, if not larger than, the quoted statistical errors. So if indeed the uncertainty due to 
calculations of the white dwarf cooling is as large as the observational error, then the LTB-ACDM model is within 
la agreement with observations. 

The age based on evolution of compact binaries is somewhat lower. In [24|[25| the authors give a result of 11.8 ± 
0.6Gyr and 11.10 ± 0.67Gyr, respectively. Moreover, the oldest known star in the Milky Way, HE 1523-0901, is 
reported to have an age of 13.2±2.7Gyr [26|, for a lower limit of 10.5Gyr. Obviously, these measurements are all in 
consistent with t^TB and ^acdm obtained with data from [l7j at la CL. Therefore, although the obtained t^TB is 
about 2Gyr younger than £acdMj it is still in la agreement with current astronomical observations, and we are not 
able to argue against the reconstructed LTB-ACDM model. 
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0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 0.36 

FIG. 1: Parameter space in Q m and Ho plane. The green and red regions stand for £acdm >11.2Gyr and £ltb > 11.2Gyr 
respectively. The black shadow region represents the la constraint to fl m and Ho from the seven-year WMAP observations 
(degeneracy is ignored). The blue region is a 2a constraint from a joint analysis from the Constitution supernovae sample, 
baryon acoustic oscillations and the five-year WMAP observations. 

We show the situation in Fig. 1. The green and red regions represent parameters with cosmic age older than 
11.2Gyr in ACDM and LTB-ACDM model, respectively. The black shadow region is a la CL constraint from the 
seven-year WMAP observations [H (we ignore degeneracy) . The blue region is a 2a CL constraint from a joint 
analysis performed in one of our previous works |27| . in which we used the Constitution supernovae sample [28| . the 
baryon acoustic oscillations [2!| and the five-year WMAP observations [3(3]. Since current limit to the cosmic age 
from astronomical measurements generally gives a result t < 11.2Gyr, we plot the regions of £ltb < 11.2Gyr (red 
shadow) and £acdm < 11.2Gyr (green shadow) in this figure. It is obvious that the ACDM perfectly passes the cosmic 
age text, while the overlap of the blue region, the black shadow region and the red shadow region implies that the 
LTB-ACDM model is also consistent with current observations. 

Finally, at the end of this section, we stress that in this paper we only consider a particular LTB model - namely, 
the one with ACDM features. Other inhomogeneous models (such as the void models) may have larger age at the 
origin, and in these cases one would have to seek for other methods to test and identify them. 

IV. SUMMARY 

In this paper we calculate the cosmic age at r = in the LTB-ACDM model, which reproduces the luminosity- 
distance and light-cone matter density of the homogeneous ACDM model. Using the constraints of O m and Ho 
from the seven-year WMAP observations combined with other cosmological observations, we get the upper limit 
tLTB < 11.7Gyr at la CL. This result is about 2Gyr younger than the cosmic age in ACDM scenario. Since current 
astronomical measurements generally put a la CL lower limit on the age of the universe of about 11.2Gyr, the LTB- 
ACDM model is still in la agreement with all these observations. However, due to the relatively younger age the 
LTB-ACDM model might be disfavored by future observations. 

Besides, even if the LTB-ACDM model successfully passes all the tests of future observations, there might be some 
other problems in this scenario. (The discussions of these complicated problems are beyond the scope of this paper.) 
The main reason is that it is difficult to fit this model into the larger framework of fundamental physics such as particle 
physics, general relativity, astrophysics, and cosmology. For example, in [l2T ] the authors mentioned the theory of 
structure formation and the integrated Sachs- Wolfe effect: study of these issues is very difficult in the LTB models. 

The topic of distinguishing the homogeneous ACDM model and the reconstructed inhomogeneous LTB-ACDM 
model is scientifically interesting and important, since it involves the question of the mysterious feature of dark 
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energy and whether the nearby region of the universe is homogeneous. This topic should be carefully investigated. In 
this paper we propose the possibility of distinguishing the LTB-ACDM scenario with the standard ACDM model by 
performing the cosmic age test. The procedure is convenient and straightforward. Although the result shows that the 
LTB-ACDM is still in la agreement with current astronomical observations, since with the same set of parameters 
this model always has a younger age than the standard ACDM model it is possible to distinguish them from future 
observations. In all, the issue of using the cosmic age test to distinguish the reconstructed inhomogeneous LTB models 
from the homogeneous ACDM model is worth further investigation, and should be taken into consideration in future 
works, e.g., in the cases when people try to construct a new model in the LTB scenario to explain the apparent cosmic 
acceleration. 
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